ON A RECIPROCITY LAW FOR FINITE MULTIPLE ZETA VALUES 



MARKUS KUBA AND HELMUT PRODINGER 

Abstract. It was shown in [7, 9] that harmonic numbers satisfy certain reciprocity 
relations, which are in particular useful for the analysis of the quickselect algorithm. The 
aim of this work is to show that a reciprocity relation from [7, 9] can be generalized to 
finite variants of multiple zeta values, involving a finite variant of the shuffle identity for 
multiple zeta values. We present the generalized reciprocity relation and furthermore a 
simple elementary proof of the shuffle identity using only partial fraction decomposition. 
We also present an extension of the reciprocity relation to weighted sums. 



1. Introduction 

Let H n = Y^k=i denote the n-th harmonic number and = Y^k=i V^ s ^he n-th 
harmonic number of order s, with n,s£N and H n = Hn^ ■ Kirschenhofer and Prodinger [7] 
analysed the variance of the number of comparisons of the quickselect algorithm [6], and 
derived a reciprocity relation for (first order) harmonic numbers. Subsequently, the reci- 
procity relation of [7] was generalized in [9], where the following identity was derived. 

j rr(a) N+l-j rr(b) 

E n N-k , \^ n N~k A i TT(b)TT(a) R (a,6) /i\ 

k=l k=l J ' 

where R^'^ = Y2k=i which can be evaluated into a finite analog of the so-called 

Euler identity for ((a)((b) stated below, 

R^ = j2( l+ b b _~ 2 )c N ^ + b-l,a+l- i ) + j2( l + a _~ 2 )Ui + a-l,b+l- i ), 
i=i ^ ' i=i ^ a ' 

(2) 



where the multiple zeta values [1, 2, 3, 5, 4, 12], and its finite counterpart are defined as 
follows: 

C(a) = C(«i, • ■ • , O := Y] ai al ~, 

n 1 >n 2 >->n r >l 1 z T 



Civ(a) = CjvK, . . . , Or) := 



n^n? 2 ...ntr 

N>m>n2>->n r >\ 1 z r 
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Note that Cn(o) = Hn ■ w = Y7i=i a « denote the weight and d = r the depth of (finite) 
multiple zeta values. The aim of this note is to derive a generalization of the reciprocity 
relation (1), stated below in Theorem 1, by considering the more general sums 



ECfc-i(&2, • • • , &a)Gv-fc(ai, • • • , a r ) (fc_i(d2, • • • , a r )(N-k(bi, . . . ,b rJ 

k=l k=t 

instead of the previously considered sums Y^k=i ~ J j^ ± an d Yl!k=i~^ 1a k ■ Our general- 
ization involves a finite variant of the shuffle identity for multiple zeta values, for which 
we give an elementary proof using partial fraction decomposition. We discuss the close 
relation between this finite variant of the shuffle identity and the shuffle identity for gener- 
alized polylogarithm functions; we will see that the finite variant of the shuffle identity is 
equivalent to the shuffle identity for generalized polylogarithm functions. To simplify the 

presentation of this work we will frequently use the shorthand notations a ( (1\ . . . . , Qi r j , 

a 2 = ( a 2, ■ ■ ■ ,a r ) and b = (b±, . . . ,b s ), b 2 = (b 2 , . . . , b s ), respectively, with r, s G N and 
ai, bk G N for 1 < i < r and 1 < k < s. 

2. The reciprocity relation for finite multiple zeta values 

We will state our main theorem below, and subsequently discuss its proof and the precise 
definition of the shuffle relation for multiple zeta values. 

Theorem 1. The multiple zeta values C7v( a ) = (N( a i, ■ ■ ■ > a r ) and Civ(b) = Civ(&i> ■ ■ ■ ,b s ) 
satisfy the following reciprocity relation. 

ECfc-i(&2, • • • , &s)Cjv-Jfc( a i> ■ ■ ■ ,a r ) y^** Cfc-i( a 2, • • • , a r )(N-k(bi, . . . , b r ) 
fc=i fe=i 

= CAr + i-,-(a)0(b) - ^ l(iV + 1 _ J j)ai + i?v(a; b). 

i^ere i?jv(a;b) = X^feLi ^~^ b ^-i (a2 '---' QV ) = R N (b;3L) satisfies a finite counterpart of the 
shuffle identity £(a)£(b) = £(aujb) for the multiple zeta value, Rn{&', b) = / jv( aLJJ b). 

Corollary 1. VFe obtain the complementary identity 

E 1 Cfc(b)Cv-fc-i(a 2 ) 5p Cfc(a)Cjv-fc-i(b 2 ) = Cj-i(b)Cv-j(a 2 ) Cv-j(a)Cj-i(b 2 ) 
{N-k) a ^ ^ (iV-A;) 6 i ~ (iV + l-j)^ J &1 

/• / \ > /i \ i Cj-i(b 2 )Cn-j(a 2 ) , . 

- Gv + w(a)0(b) + ^ i(iy + 1 _ J j)ai + ^v(a; b). 

Next we state an immediate asymptotic implication of our previous result. 

Corollary 2. For N — 2n + 1, j — n + 1, with a±, b\ G N \ {1} and n — > oo we obtain the 
following result. 

Um / ^ + ' C fc -i(a 2 )Cv- fc (b) \ = 2C(a)c(b) . 

v fc=l fc=l 7 
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In order to prove Theorem 1 we proceed as follows. 

Cfc-i(t>2)0v-fc(a) r ( v Cfc-i(b 2 ) ^ O-i (a 2 ) 

fc=i fc=i e=N+i-j 

N-l t , s N-l . , u x 

= Civ-i(a)C-(b) + 2_ m 2- 



J=JV+l-j fc=l 



/■ fU\ i 0~i(a2)Gv-*(b) 
Cjv-i(a)0(b) + ^_ ^ 

£=N+l-j 

JV 

CiV+w(a)Q(b) + ______ + 2^ 

=7V+2-j 



Cjy-j(a 2 )Cj-i(b 2 ) 0-i(a 2 )Qv-*(b) 
(_/V + 1 — j) ai j bl ^— ' 



This proves the first part of Theorem 1 and 



D / ,n CA r -fc(b)Cfe-i(a2, • • • , Or) 
«iv(a; b) = 2^ 7^ • 



k=l 



For the evaluation of Rn{&] b) we note that -Ro( a ; b) = 0, and further 



A' 



R N (&; b) = J2 {M*, b) - i4-i(a; b)). (3) 



fc=i 

We have 

R N (a; b) - i^a; b) = ^ _____ 

Now we use partial fraction decomposition 1 , which appears already in [11], 



a (i+b-2\ b H+a-2\ 

V 1 b ~ l ' + V 1 a ~ l ' (A) 

L^t Afi+b-lUa+l-i Mi+a-1 ( AT _ U\b+l-i ' V ) 



k a {N - k) b ^ N i +b -i k a+i-i Z_, N i+a-i( N _ k y+i 
and obtain 

E 1 Cfc-l(>2, «rRjV-l-fc(fr2, _A) _ "l 2 ) Cfc— 1 («- , • • • , ar)Qv-l-fc(&2, ■ ■ ■ A) 

C_/Y _ fc)bifcai _^ ]yi+6i-lj^oi+l-j 

fc=l V y i=l fe=l 

^ ^ (^-l 2 )^" 1 ^ 2 ' ' ' ' ' a r)0v-l-fc(&2, ■■■,b s ) 



jyi+ai-lfjy _ nSi+l-i 
i=l fc=l v ; 



This identity has been rediscovered many times. For a fascinating historic account, see [8]. 
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Consequently, by summing up according to (3) we get the following recurrence relation for 
i?Ar(a; b). 

R N (a; b) = ^ ^2 i+bi-i ^ni-i(»i + I - i,a 2 , . . . ,a r ; b 2 , 

i=l ni=l 1 

+ ^2 i+ai-l ^i-l(°2» • • • , a r\ h + 1 - 1, b 2 , 

i=l ni=l 1 

This recurrence relation suggests that there exists an evaluation of -Rat (a; b) into sums of 
finite multiple zeta values, all of them having weight w = Y7i=i a r + YH=i b% an d depth 
d = r + s. In order to specify this evaluation we need to introduce the shuffle algebra for 
(finite) multiple zeta values. 

2.1. The shuffle algebra. Let A denote a finite non-commutative alphabet consisting of 
a set of letters. A word w on the alphabet A consists of a sequence of letters from A. 
Let A* denote the set of all words on the alphabet A. A polynomial on A over Q is a 
rational linear combination of words on A. The set of all such polynomials is denoted by 
Q{A). Let the shuffle product lu be defined on < Q{A) as follows: for any w, v G A* with 
w = X\ . . . x n , v = x n+ i . . . x n+m , Xi G A for 1 < i < n + m 

WluV := ^ X a(l)Xa(2) ■ ■ ■ X<r(n+m) , (6) 

where the sum runs over all (™^ m ) permutations cr G & n +m which satisfy o - " 1 ^) < O"" 1 ^) 
for all 1 < j < k < n and n + l<j<k<n + m. The sum is over all words of length 
n + m, counting multiplicities, in which the relative orders of the letters X\, . . . ,x n and 
x n+ i, . . . ,x n+m are preserved. The term "shuffle" is used because such permutations arise 
in riffle shuffling a deck of n + m cards cut into one pile of n cards and a second pile of m 
cards [4]. Equivalently, we can recursively define the shuffle product as follows. 

Vw G A*, €luw = WLue = w, 

Vx, y G A, w, v G A*, xwluijv = x(wlu?/v) + i/(iwluv). (7) 

2.2. The shuffle algebra and multiple zeta values. Let a and b denote the multi- 
indices a = (a 1; . . . , a r ) and b = (bi, . . . , b s ) with dj, bj G N for 1 < i < r, 1 < j < 
s. To any multi- index we associate a unique word over the non commutative alphabet 
A = {u; ,^i}. Let A = A(a), B = B(b) such that A := uj^ujxuj^ujx . ..cu^^ui and 
B := Uq^ujiuJ^^uji . . .lu^^lui. To each word we associate a finite multiple zeta values 
by the following linear correspondence: Ztv^i^q" 1 ) = C7v( a )> an d in general 

lu^ 1 ^) = Cjv(oi, • • - ,a r ). 

(8) 



b„ 



(5) 



Z( a\ — 1 ai — 1 a r — 1 \ \ v ry i a.2 — 1 

n[u l loiUq u)i...u Q r lox J = ) j — Z rai _i(o; ^ u x ... 

N / 72/1 

m=l 1 
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Moreover, assuming that = (a^i, . . . , at, re ) and with 1 < I < h, qt e Q, and Ae = A(a^), 
we get by linearity of the correspondence 

(h \ h h 

^2 Qe A e ) = E qeZ N {A e ) = ^ ftCivfc, • • • , a i>n ). 
i=i ' e=i £=i 

We observe that the partial fraction decomposition (4) of ka ^_ k ^b above mimics the 
basic shuffle identity for words A = Uq~ 1 lj 1 , B = uj b ~ l uji, 

A^B = |j ( 6 ~ l _ | ty^u^- 1 -^! + g ^ ~- 1 ^o- 1 ^^ 6 - 1 -^, 

which appeared in Hoang and Petitot [10]. 

The key to the explicit evaluation of -Rjv(a; b) is the following result concerning the 
shuffling of the words A and B, associated to the multi-indices a and b. 

Lemma 1. Let A := uJq 1 ~ 1 ujiuJq 2 ~ 1 uji . . . ujQ r ~ 1 uJi and B := uJq 1 ~ 1 ujiuJq 2 ~~ 1 uji . . . uJq s ~ 1 uji, 
with a = (ai, . . . , a r ) and (b\, . . . , b s ) with a^, bj£N,l<i<r,l<j<s. We have 

A - B = E f + b -7 2 W bi ~V(^5 2 ) + E f + ai _~ 2 V--v(a 2 ^o, 

i=l ^ ^ ' i=l ^ ^ ' 

with A\ := uj^ujiujI^uji . . . lu^' 1 ^, B[ := ^"Vc^ 2 " 1 ^ . . . u^ui and further A 2 : = 
Uq 2 ^ . ..lo^ujx, B 2 := . . .^ S ~V- 

The special case r = s = 1 is a result of Hoang and Petitot [10]; we simply use the 
recursive definition of the shuffle product (7) and obtain the result of Lemma 1. 

Now we are ready to prove the evaluation of -Rjv(a;b). Let A = A(a) and B = A(h) 
denote the words associated to the multi-indices a and b, 

R n (sl; b) = Z N (A^B) =: Ov(amb). (9) 

We use induction with respect to the depth d = r + s. The result clearly holds for d = 2, 
(2), as shown in [9]. Assuming the result for all depth r + s < d we obtain according to 
the induction hypothesis 

R N (&; b) = V" i+^-i flm-i( Q i + 1 - h «2, • • • , a r ; b 2 , . . . , b s ) 

i=i m=i 1 

+ E E i+ai-1 ^m-l( Q 2, • • • , Or! 6l + 1 - 1, fa, ■ ■ ■ , b s ) 

t=l n±=l 1 

= E E i+bi-l ^m-lCALuSa) + ^ E ^i+ai-l ^m-l(^2LU-Bi). 



1 1 ^1 1 1 ^1 

i=l m=l 1 i=l ni=l 1 
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By (8) and Lemma 1, using the notations of Lemma 1, we get 

Zn(A^B) = ^ i+fti-l Z m-l(A^B 2 ) +YY1 J+ai-1 Z ni -i(A 2 ^Bi). 

1=1 ni=l 1 i=l ni=l 1 

Consequently, 

#jv(a; b) = Z n (AhaB) = Ov(a^b). (10) 
This proves the stated result for -Rat (a; b) and the remaining part of Theorem 1. Moreover 
we observe that the depths d = r + s and the weights w = Yll=i a « + J2t=i of the finite 
multiple zeta values are all the same. Corollary 1 can easily be deduced by noting that the 
sum of the left hand sides of Corollary 1 and Theorem 1 add up to -R n (a; b) with respect 
to two extra terms. Now we turn to the proof of Corollary 2. Since for N = 2n + 1 and 
j = n + 1 and n->oowe have 

lim Cj(b)CN+i-j(b) = lim C n+1 (b)C„+i(b) = C(a)C(b), 

( n (b 2 ,...,b s )C n (a 2 ,...,a r ) 
lim — — -7 = 0, 

(n + 



n— >oo 



lim R N (&;b) = lim C2n+i(aLub) = C(aujb) = C(a)C(b), 



n— +oo 



the result of Corollary 2 can be immediately deduced from Theorem 1. Note that the last 
identity is the well known shuffle identity for multiple zeta values; we refer the reader to 
the excellent article [2]. 

3. POLYLOGARITHMS AND THE FINITE SHUFFLE IDENTITY 

Let Li a (^) = Li olj .„ )0r (z) denote the (multiple) polylogarithm function with parameters 
Cb\ , . . . , CL r , defined by 



Lia(z) =Li aii ... iap (z) = V ^nr 2 -, (11) 

z — ' Ut fin . . . n°: r 

ni>n 2 >-->n r >l 1 z r 

We note that the value -Rat (a; b) can be obtained in the following way. 

D / ,n Cfc-i( a 2, • • • > a r )(N-k{bi, ■ ■ ■ , b s ) r Ar-,Li a (z) Li b (^) 
R N {&; b) = 2^ ^ = [z \ YZ~ Z 

fe=i 

Consequently, our finite shuffle identity (9) for i?^( a ib) is equivalent to the following 
shuffle identity for polylogarithm functions. 

Lia(z) Li b (z) = U aLUh (z). 

Note that by evaluating at z — 1 the shuffle identity for polylogarithm functions implies the 
shuffle identity for multiple zeta values. The identity above is well known, see for example 
the article [2]. The shuffle identity for polylogarithm functions is due to the iterated 
Drinfeld integral representation of polylogarithm functions and multiple zeta values due 
to Kontsevich [12]. As remarked in [2] the shuffle product holds since the product of two 
simplex integrals consists of a sum of simplex integrals over all possible interlacings of 
the respective variables of integration. Our result for -Rat (a; b) implies that the shuffle for 
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polylogarithm functions, and also for multiple zeta values, can be derived using only basic 
partial fraction decomposition. 

4. The reciprocity relation for weighted multiple zeta values 

Results similar to Theorem 1 and Corollary 2 can be obtained for products of weighted 
finite multiple zeta values, Cn(cli, a 2 , . . . , a r ; a±, . . . , ay), a, G R \ {0} for 1 < i < r, defined 
as follows: 

Civ(a, cr) = Cjv(oi, a 2 , ■ . ■ , a r ; a x , . . . , ay) — ^ z= — . 

JV>m>n 2 >-'>n r >l Ali=l i i 

Of particular interest are the cases 0{ G {±1} corresponding to a mixture of alternating 
and non- alternating signs. We only state the result generalizing Theorem 1, with respect 
to the notations a 2 = (a 2 , . . . , a r ), er 2 = (a 2 , • • • , ay) and the corresponding notations for 
b 2 and r 2 , and leave the generalizations of Corollaries 1,2 to the reader. 

Theorem 2. The multiple zeta values Gv( a 5 CT ) o,nd Civ(k, T ) with weights a and r satisfy 
the following reciprocity relation. 

E Cfc~i(b 2 ,T 2 )Cjv_ fc (a,o-) N ^ 3 Cfc-i(a 2 ,o- 2 )Gv_ fc (b,T) 

k=l 1 k=l 1 

> / \ ^ / 1 \ 0-i(b2, T 2 )CAr-j(a 2 , cr 2 ) 

r iJ °~i (i\+l—j) ai 

Here i?iv(a, cr; b, r) = Xlfc=i ^~ fc ^ b '^l^ L 1 ^ a2 ' q ' 2 ' > = -Rjv(b, r; a, cr) satisfies an analogue of 
the shuffle identity with respect to the weights cr and r. 

The proof of Theorem 1 can easily be adapted to the weighted case. Hence, we only 
elaborate on the main new difficulty, namely the evaluation of the quantity 

o / ,N Ov-fc(b, r)Cfc-i(a 2 ; cr 2 ) 
i2jv(a, cr; b, r) = ^ -j— . 

k=l 1 

Proceeding as before, i.e. taking differences and using partial fraction decomposition, we 
obtain the recurrence relation 

ai N ^i+6i-2\ 

R N (&, cr; b, r) = Y] Y] i+llZl — -Rm-iCfli + 1 - i, a 2 , — , cr 2 ; b 2 , r 2 ) 

i=l ni=l 1 1 
6i AT 

+ .i+al-U ^'- 1 ^ 2 ' CT2; 6l + 1 ~ b2 ' ^' T2) - 

i=l ni=l 1 1 

Consequently, the value (a, cr; b, r) can be evaluated into sums of weighted finite mul- 
tiple zeta values according to a shuffle identity with respect to the weights cr and r. We 
omit the precise definition of this generalization and leave the details to the interested 
reader. 
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Conclusion 

We presented a reciprocity relation for finite multiple zeta values, extending the previ- 
ous results of [7, 9]. The reciprocity relation involves a shuffle product identity for (finite) 
multiple zeta values, for which we gave a simple proof using only partial fraction decom- 
position. Moreover, we also presented the reciprocity relation for weighted finite multiple 
zeta values. 
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